Huygens-Fresnel-Kirchhoff construction for quantum propagators with application to 

diffraction in space and time 
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We address the phenomenon of diffraction of non-relativistic matter waves on openings in absorb- 
ing screens. To this end, we expand the full quantum propagator, connecting two points on the 
opposite sides of the screen, in terms of the free particle propagator and spatio-temporal properties 
of the opening. Our construction, based on the Huygens-Fresnel principle, describes the quantum 
phenomena of diffraction in space and diffraction in time, as well as the interplay between the 
two. We illustrate the method by calculating diffraction patterns for localized wave packets passing 
through various time-dependent openings in one and two spatial dimensions. 

PACS numbers: 03.65.Nk, 03.75.-b, 42.25.Fx 



I. INTRODUCTION 

Diffraction and interference of matter are among the 
most fascinating and controversial aspects of quantum 
theory. It is not surprising that laboratory exploration 
of these phenomena with subatomic, atomic, and molec- 
ular particles has been at the heart of experimental re- 
search since early days of quantum mechanics [l| . As of 
today, a wave-like behavior of matter has been success- 
fully demonstrated for a number of elementary particles, 
atoms, simple molecules (see Ref. [l| for a comprehen- 
sive review), and, most notably, for some heavy organic 
compounds including Ceo 0, C70 [H, and C60F48 0|. A 
majority of these experiments involve sending a mono- 
energetic beam of particles through a screen with open- 
ings (apertures) such as slits, diffraction gratings, or Fres- 
nel zone plates. Mathematically, the problem of quan- 
tum diffraction on stationary spatial apertures can be 
described by the Poisson equation and treated by meth- 
ods originally developed in the context of diffraction and 
interference of light 0, @ . 

A diffraction phenomenon of a different kind - "diffrac- 
tion in time" - was introduced by Moshinsky six decades 
ago and subsequently studied by many researchers, 
both experimentally § and theoretically (see Refs. (9l.[l0j 
for reviews) . The phenomenon has to do with time evolu- 
tion (more precisely, with non-uniform spreading) of ini- 
tially sharp wave fronts in quantum systems, and mani- 
fests itself already in one dimension. Thus, in the original 
set-up proposed by Moshinsky Q, a perfectly absorbing 
shutter is placed in the way of a mono-energetic beam of 
non-relativistic quantum particles. Then, a sudden re- 
moval of the shutter creates a "chopped" particle beam 
with a sharp wave front. As shown by Moshinsky, such a 
wave front disperses non-uniformly in the course of time 
and, most interestingly, develops a sequence of diffrac- 
tion fringes. Mathematically, these fringes appear to be 
analogous to the ones observed in diffraction of light on 
the edge of a semi-infinite plane. 

There are several ways of treating quantum diffrac- 
tion theoretically. One commonly used, physically moti- 



vated method for evaluating the wave function of a quan- 
tum particle passing through an opening in a diffraction 
screen is the "truncation" approximation, which is based 
on the composition property of quantum propagators. In 
this approximation, transmission of a spatially localized 
wave packet through the diffraction screen is treated as 
a three stage process: (i) the wave packet is propagated 
freely during the time that it takes the corresponding 
classical particle to reach the screen, then (ii) the wave 
function is reshaped (or truncated) in accordance with 
the geometry of the aperture, and, finally, (iii) the resul- 
tant wave function is propagated freely for the remaining 
time interval. The reader is referred to Refs. (HUH] for 
details and implementation examples of the truncation 
approximation. 

One drawback of the truncation approximation is that 
it does not account for diffraction in time. Brukner and 
Zeilinger fl5j proposed another, more versatile method 
for solving the problem of quantum diffraction. The 
central assumption of their method is that the time- 
dependent wave function in question satisfies certain ex- 
plicitly known, inhomogeneous time-dependent Dirich- 
let boundary conditions at the surface of the diffraction 
screen. More specifically, two assumptions are made: (i) 
the value of the wave function at a point inside the aper- 
ture is assumed to equal the value that the wave func- 
tion would have at this point if the diffraction screen 
was absent, and (ii) the wave function vanishes at every 
point of the diffraction screen outside the aperture. The 
first assumption is commonly referred to as the Kirch- 
hoff approximation, while the second corresponds to the 
physical assumption of perfect reflectivity of the screen. 
The method of Brukner and Zeilinger has since been suc- 
cessfully used by several researchers to study quantum 
diffraction in both space and time (l6l - [l8j . 

In this paper, we present another approach to quan- 
tum diffraction in space and time. Our method is based 
on the Huygens-Fresnel principle and Kirchhoff theory 
of diffraction, and allows one to calculate the time- 
dependent quantum propagator for the problem of parti- 
cle diffraction on spatio-temporal openings in otherwise 
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perfectly absorbing screens. Our expression for the prop- 
agator is especially adaptable for calculating diffraction 
patterns in situations in which the initial wave function 
is given by a spatially localized wave packet or by a su- 
perposition of several localized wave packets. Similar to 
the method of Brukner and Zeilinger, our construction 
provides a unified framework that treats the phenom- 
ena of diffraction in space and diffraction in time on the 
same footage. The main difference between the method 
of Brukner and Zeilinger and our method is that the 
former is designed to treat perfectly reflecting screens, 
while the latter is only applicable to perfectly absorb- 
ing ones. Finally, our approach allows for calculation of 
quantum diffraction patterns produced by openings in 
spatially curved screens. 

The paper is organized as follows. In Sec. [TT| we 
develop an expansion of the propagator for a non- 
relativistic quantum particle passing through a time- 
dependent opening in an absorbing screen. In Sec. IIIII 
we demonstrate utility of the expansion by applying it to 
some example systems in one and two spatial dimensions. 
In Sec. II VI we discuss our results and make concluding re- 
marks. Some technicalities are deferred to an Appendix. 

II. HUYGENS-FRESNEL-KIRCHHOFF 
CONSTRUCTION FOR QUANTUM 
PROPAGATORS 

In this section, we address two quantum-mechanical 
phenomena - diffraction in time, pioneered by Moshinsky 
0, H, 03 1 and diffraction in space, as described by Kirch- 
hoff theory 0, @. We recast standard descriptions of 
both phenomena in a way analogous to a time-dependent 
formulation of the Huygens-Fresnel principle. Expressed 
in this way, the two diffraction processes appear to be 
closely related and can be straightforwardly unified into 
a single model of diffraction in space and time. 

A. Notation 

In order to facilitate the clarity of the following presen- 
tation, we begin by introducing central physical quanti- 
ties and fixing notation. 

The focus of this paper is on the quantum propagator 
K(q, q';t) that describes the motion of a quantum par- 
ticle in the /-dimensional coordinate space by relating a 
particle's wave function \&(q;t) at time t to that at time 
t = through 

*(<«<)= / dV^(q,q';*)*(q';0). (1) 

JRf 

The propagator is the solution of the time-dependent 
Schrodinger equation 

dK 

ih— = HqK (2) 



with the initial condition 

lim#(q,q';t) = 5(q-q'), (3) 

where Hq denotes the Hamilton operator in the position 
representation. (The propagator is also subject to certain 
absorbing boundary conditions that we do not discuss at 
this point.) In the case of a free particle of mass m we 
have Hq = ~^y\ and K = X (q - q'; t) with 

Hereinafter, the subscript V indicates that the corre- 
sponding quantity refers to the case of a free particle. 

We also consider the energy-domain Green function 
defined as the Laplace transform of the propagator, 

/>oo 

G(q,q';£) = / die-*Jir(q,q';<) 
Jo 

= £[X](q,q'; S ) with s = ^. (5) 

in 

Consequently, the propagator is obtained from the Green 
function by means of the inverse Laplace transform, 
K = £ _1 [G]. In the free particle case, Eqs. (U]) and (O 
define the free-particle Green function Go(q — q';E) = 
C[K n ](q — q'; s) that satisfies 

v*Go+fc 2 G = -f?<Kq-q') with fc 2 = ^.(6) 

B. Diffraction in time 

We now address the phenomenon of diffraction in time, 
first considered by Moshinsky Ml a nd later explored by 
many researchers (see Refs. [9|, LL(| for reviews). In its 
one-dimensional formulation, the Moshinsky problem is 
concerned with time evolution of a quantum particle, 
whose wave function ^(£;t) is localized to the semi- 
infinite interval (— oo,a;i) at time t = 0, i.e., ^(^0) = 
for £ > x\. Over time, an absorbing wall (shutter) 
is switched "on" and "off" at the point x\, according 
to a protocol defined by a characteristic function x(t)- 
The latter is allowed to take values between and 1, 
with representing the case of perfect absorption (shut- 
ter "on") and I corresponding to perfect transmission 
(shutter "off"). The open interval < \ < 1 repre- 
sents the case of partially absorbing and partially trans- 
mitting (but reflection-free) shutter. One is then inter- 
ested in the wave function ^>(^;t) of the particle to the 
right of the shutter, £ > x\, at t > 0. Moshinsky Q 
analyzed this problem for a "monochromatic" incident 
wave \&(£;0) = 6(2:1 — £)e 1 ^, where k > and 9 is 
the Heaviside step function, and a perfectly absorbing 
shutter that gets suddenly removed at an instant to, i.e., 
x(t) — 0(i— to). His analysis showed that at times t > to, 
the front of the probability density wave exhibits patterns 
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identical to those observed in the Fresnel diffraction of 
light from the edge of a semi-infinite plane, thus giving 
rise to the term "diffraction in time" . 

Our objective is to devise an expression for the propa- 
gator K(x,x';t) that describes Moshinsky diffraction by 
a shutter positioned at a point x\ , such that 



x < X\ < x , 



(7) 



and controlled (opened and closed) in accordance with 
the characteristic function x(t) of an arbitrary functional 
form. Our construction relies on the Huygens-Fresnel 
principle, which, for the purpose of the current problem, 
can be formulated as follows: The disturbance at the 
point x produced by a source O' , located at some other 
point x', can be viewed as produced by a fictitious source 
0\, located at a point x\ in between x' and x, cf. Eq. (J7J. 
The strength of the fictitious source 0\ is determined by 
the disturbance at the point X\ produced by the original 
source O'. When a (partially) absorbing shutter is placed 
at the point xx, the strength of the fictitious source 0\ 
is modulated by the characteristic function x(t) taking 
values between and 1. Mathematically, this can be 
summarized as 

K(x,x';t)= / dti u(x — xi, x\ — x'\ t, ti) 
Jo 

x K (x - xr,t - ti)x(h)K (xi -x';ti). (8) 

where u, having dimensions of speed, is a yet-to-be- 
determined function of the distances x — x\ and x\ — x' 
and times t and t\. 

The physical meaning of Eq. © is transparent: The 
probability amplitude of an event in which the particle 
goes from x' to x in time i can be expressed as a sum of 
probability amplitudes over all composite events in which 
the particle first goes from x' to an intermediate point x\ 
in a time t\ and then reaches x from x\ in the remaining 
time t — t\. Because of their simple physical interpreta- 
tion, propagator expansions similar to Eq. (JSJ are often 
used for qualitative description of quantum interference 
phenomena (e.g., see @ for a qualitative discussion of 
a two-slit interference experiment), however the explicit 
functional form of u is usually neither specified nor taken 
into account. 



It is interesting to observe that Eq. (j9)) does not specify 
the function u uniquely. In fact, Eq. @ turns out to be 
an identity satisfied exactly by infinitely many different 
functions u, some examples being (see App. [A")) 



X X\ , , X ^ X 
U = J]- — + (1 - 1]) 7 



t - u 



2ih exp(— ( 2 
nmt erfc(C) 



. , >9 mix — x') 2 

with C = — i '— 

2iht 



(10) 
(11) 



Here, ry is an arbitrary complex number, and "erfc" de- 
notes the complementary error function. 



2. Moshinsky shutter, x(t) = 0(t — to) 

We now show that the functional form of u can be 
uniquely determined by comparing Eq. ((5J) with an ex- 
act expression for the quantum propagator in the original 
Moshinsky set-up 0], in which the absorbing shutter, lo- 
cated at xi, is closed until a time to and open afterwards, 
i.e., x(*i) = 9(*i - to) for < t , ti < t. 

On one hand, a direct construction of the propagator 
Ky[(x,x';t) for the original Moshinsky problem leads to 



K M (x,x':t) = J dx" K (x-x";t-to)K (x" -x';t ) 
= Kq(x — x ; t) 



1 — Trerfc [ (X\ — In) ( / —- ; r 

2 l v u; V 2iht (t-t ) 



with 



to . it — to 

Xq = x h x 

t t 



(12) 



(13) 



The first equality in Eq. (|12[) combines the composition 
property of quantum propagators and the fact that, at 
time to, all the probability density to the right of x\ has 
been absorbed by the shutter. This leads to the trunca- 
tion of the upper limit in the x" integral. 

On the other hand, a substitution of x(ti) = 0(ti — to) 
into Eq. © yields 



1. Free particle, x(t) = 1 

Interested in determining the functional form of u = 
u(x — x\,xi — x';t,ti), we first direct our attention to 
the simplest possible scenario, in which the shutter stays 
open throughout the entire time interval from to t, i.e., 
x(ti) = 1 for < t\ < t. In this case, the propagator in 
the left-hand side of Eq. © is given by the free particle 
propagator, K = K , requiring u to satisfy 

Kq(x — x'; t) = / dti u Kq(x — xi;t — ti)Ko(x\ — x' ;ti) . 
Jo 

(9) 



Kyi(x,x';t)= I dti u(x — xi, x\ — x'; t, ti) 



x K {x-xx;t-tx)K Q (xx -x';ti). (14) 



Equations (fT2j) and (fT4j) allow us to uniquely determine 
the function form of u by requiring Km = Km- Indeed, 
let us for the moment fix the values of x, x' , and t, and 
treat the propagator Km as a function of the shutter 
opening time to only, i.e., Km = ^M(to)- First, we note 
that lim to _,. 4 - Km = 0. Indeed, xq — > x as to ~> t~, 
and the argument of the complementary error function 
in Eq. (jT2"j) tends to e i37r / 4 oo, making the value of the 
complementary error function approach 2. This limit, of 
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course, corresponds to the trivial case of the absorbing 
shutter being closed throughout the entire time interval 
from to t. Then, in view of this limit, we rewrite the 
Moshinsky propagator as 

£ M (to) = -/V^%^. (15) 

A straightforward (but somewhat tedious) calculation 
yields 

di^M^i) 1 fx — X\ X\ — x'\ 
ddl ~ ~ 2 V t - h H h ) 

x K {x - x x ;t - ti)K (xi -x';ti). (16) 

It is now clear that the Moshinsky propagator Km, as ex- 
pressed by Eqs. ([15)1 and (fll)|). is equivalent to the prop- 
agator Km, given by Eq. ([14")) with 

1 fx- Xi Xi - x'\ 

u= 2{j^ + —)' (17) 

Note that Eq. (JT7]| is equivalent to Eq. (JIDJ with rj = 
1/2. Also, Eq. ([i"7]) provides the physical meaning of the 
function u: The latter is a characteristic (mean) velocity 
of the particle when it traverses the shutter. 

3. Arbitrary x(i) 
A substitution of Eq. dTTJ) into Eq. (gj) yields 

x Koix-x^t-t^xit^Koixx-x'ih). (18) 

Here, the spatial points x, x' , and X\ are subject to the 
condition given by Eq. ([7]), and the characteristic func- 
tion x is allowed to take values between (perfect ab- 
sorption) and 1 (perfect transmission). 

Equation ([T8|) provides a formulation of the Huygens- 
Fresnel principle for a one-dimensional particle in the 
presence of a point-like absorbing obstacle, whose ab- 
sorbing properties change in the course of time. At this 
point, it is important to emphasize that the analysis pro- 
vided in this section should not be regarded as a rigor- 
ous mathematical proof of the propagator expansion (fT"8)) . 
Unavoidable difficulties in solving the problem from the 
first principles stem from the lack of a proper unambigu- 
ous definition of point-like, generally partial and time- 
dependent, absorption. In our approach, however, we by- 
pass this issue by modeling the absorption with the help 
of a time-dependent characteristic function, x(t), and re- 
lying on the validity of the Huygens-Fresnel construction 
in its most general form, Eq. (JSj) . Consequently, we re- 
move any arbitrariness in the Huygens-Fresnel construc- 
tion by determining the function it, originally unknown in 
Eq. ([5]), through analyzing the case of x(t) corresponding 



to the Moshinsky shutter problem, for which point-like 
absorption can be defined unambiguously. 

It is interesting to note a formal similarity between 
Eq. (1181) and the well-known Lippmann-Schwinger equa- 
tion |9l.llll. Il9| , which in the case of a point-like perturba- 
tion, situated at x\ and described by a spatio-temporal 
potential of the form V(£,t) = d(£ — xi)U(t), reads 

K sc (x,x';t) = K Q (x - x';t) 

i r l 

--J dhKoix-x^t-t^UitJK^ixux'ih). (19) 

As before, K is the (free-particle) propagator in the ab- 
sence of the perturbation potential, and K sc is the propa- 
gator corresponding to the full scattering problem. Equa- 
tion (1191) gives rise to a multiple collision representation 
of the scattering propagator, known as the Dyson series 

Despite their superficial resemblance, Eqs. (fT5j) and 
(|T9| describe very different physical processes: The 
Lippmann-Schwinger equation pertains to the phe- 
nomenon of quantum scattering, whereas the Huygens- 
Fresnel propagator expansion represents quantum mo- 
tion in the presence of obstacles that (partially) absorb 
matter waves without deflecting them. In fact, it is easy 
to show that the Lippmann-Schwinger equation can not 
be used to model perfect absorption, corresponding to 
the trivial choice x = in Eq. (Tl8)) . Indeed, the case 
of a perfectly absorbing obstacle at xi would require 
K sc (^,x';t) = Q(x 1 — {;)K (l; — x';t), which is clearly 
incompatible with Eq. (|19[) . 

In order to avoid possible confusion, we emphasize 
that it is the assumption of a point-like perturbation, 
V(£,t) = (5(£ — xi)U(t), that does not allow the 
Lippmann-Schwinger equation to properly capture the 
physics of absorption. On the opposite, the Lippmann- 
Schwinger equation with a smooth, complex-valued po- 
tential function V(£, r), defined over an extended spa- 
tial interval, is often the method of choice in modeling 
absorbing boundaries (see Ref. [2(| for a comprehensive 
review) . 



4- Continuously opening shutter, x(t) = e 

In order to demonstrate the usefulness of the Huygens- 
Fresnel formulation of diffraction in time we apply the 
propagator expansion (TT8"]) to a modified Moshinsky shut- 
ter problem, in which the absorbing shutter is initially 
closed, x{ti) — for t\ < 0, and then opens continuously 
in accordance with 

X(h) = exp(-r/t x ) for t x > . (20) 

Here, r > determines the rate at which the shutter 
opens. Note that the shutter described by Eq. ([20]) is 
completely removed only in the limit of infinitely long 
times, lim x{ti) = 1- 

t\ — >+oo 
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Using the Huygens-Fresnel expansion (|18|) , we now de- 
rive an explicit expression for the propagator K(x,x';t) 
connecting a point x' to the left of the shutter (x' < x±) 
at time with another point x to the right of the shutter 
(xi < x) at time t. To this end, we rewrite Eq. (fl"8|). in 
view of Eq. (j4|) , as 



K(x,x';t) 



m f d d 
2m \ dx dx' 



dtiK (x - xi;t - ti)x(ii)#o(a:i - x';t x ) . (21) 



Then, taking into account Eq. (|20|). we write 
x{t\)Ko(xi — x';t\) = Kq(xi — x;ii), where x satisfies 



(xi 



x? 



(xi 



„'\2 



2ihr 



(22) 



t > 0. Second, it is clear from Eq. (|28|) that the propaga- 
tor K(x,x';t) is not symmetric in the coordinates x and 
x' . The absence of such symmetry is typical for quantum 
motion in the presence of time-dependent obstacles (e.g., 
see Sec. A.l in Ref. [1]). 



C. Diffraction in space 

We now address Kirchhoff theory of diffraction in space 
In particular, we rewrite the Kirchhoff 's formula- 
tion in its time-dependent form, in which it can be di- 
rectly applied to diffraction of quantum wave packets. 

To this end, we consider a smooth (/ — l)-dimensional 
surface S that is defined as the zero set of a real-valued 
function s : W — > R, 



Solving Eq. (|22j) for x and choosing the solution that 
converges to x' in the limit r —¥ + , we have 



x = X\ 



where 



{xi 



= — tan 
2 



+ /2|>r 

\ TO 

2Ht 



(23) 

(24) 
(25) 



to(xi — x') 2 
Equation (|21[) can now be rewritten as 

2f( '•£) f 9 dx d 

2m \ dx dx 1 dx 

f dhKoix-x^t-t^Koixi-x;^). (26) 
Jo 

The integral over t\ can now be evaluated explicitly: 

/ dt 1 K [x-x 1 ;t-t 1 )Ko(x 1 -x;ti) 
Jo 



2ih 



erfc 



TO _ 

—-—(a; — x) 

2iht v ; 



(27) 



This identity is derived in Appendix [A] for the case of 
real-valued x (see Eqs. (IA2|) and (|A7[)): the derivation 
however can be extended to the case of Rex < x\ < x 
and Imx < (cf. Eqs. (123H25jl ). Finally, substituting 
Eq. (|27l) . together with dx/dx' = (x\ — x')/{x\ — x), into 
Eq. (|26|) and taking the partial derivatives with respect 
to x and x, we arrive at 



i , 1 / Xl — x' 
2 V xi — x 



ifo(x — x; i) 



(28) 



Two quick remarks are in order. First, Eq. (|28[) guar- 
anties that, as expected, K(x,x';t) —t K Q (x — x';t) as 
r — > + . This limit corresponds to the standard Moshin- 
sky set-up in which the shutter stays completely open at 



S = { qi € : s(qi) = 0} . 



(29) 



The surface S is assumed to partition the position space 
into two disjoint regions, such that the function s takes 
different (and constant) signs in the two regions. We 
also consider two spatial points, q and q', that lie on 
the opposite sides of the surface. For concreteness we 
take s(q) > and s(q') < 0. The surface S gives 
the location of a non-transparent, absorbing screen, in 
which some transparent openings (apertures) may be 
"cut out". These openings and, more generally, the ab- 
sorbing properties of the screen, can be described by a 
spatially-dependent characteristic function x( ( li) taking 
values between zero (perfect absorption) and one (per- 
fect transmission) at a point qi e S. Kirchhoff theory of 
diffraction [5j, IfJ allows one to express the Green function 
G(q, q'; E), connecting the points q and q' at an energy 
E, as an integral along the screen: 



G(q,q';£;)=-— / d^ <5( S (q 1 )) X (q 1 ) 
2m J Rf 

x Vs(qi) • ( G (q - qi; £)V qi G (qi - q'; E) 



Go(qi-q';£)V qi G„(q-qi;£) 



(30) 



It is important to emphasize that Kirchhoff method, 
Eq. (j3"0)l . is not applicable to diffraction on apertures 
in transmission screens with reflecting (e.g., Dirichlet or 
Neumann) boundary conditions. Instead, Kirchhoff the- 
ory is known to be a good model for diffraction on per- 
fectly absorbing, or "black" , infinitely thin obstacles j2l1 — 
[23j . In fact, Eq. (|30|) provides an exact solution to the 
time-independent diffraction problem with the boundary 
conditions on the screen given by the so-called Kottler 
discontinuity (see Refs. [22|, HH and references within): 
The probability amplitude field has a discontinuity across 
the screen equal to minus the value of the free-space field 
at that point. A similar condition is imposed on the nor- 
mal derivative of the field. 
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For our purposes, it is important to construct a time- 
dependent formulation of Kirchhoff diffraction. This is 
straightforwardly achieved by first rewriting Eq. (|30|) as 

GM;E) = ~ f d/q 1 5(s(q 1 ))x(q 1 ) 

x Vs(qi) ' (V q - V,0 G (q - qi; S)G (qi - q'; E) , 

(31) 

and then performing the inverse Laplace transform from 
energy-dependent Green functions to time-dependent 
propagators, yielding 

K(q,q';t) = l f dh f dAfc $(*(qi)) 

Z JO JRf 

( q - qi , qi - q' \ „ , s 

x #o(q-qi;<-ii)x(qi)#o(qi-q';ti). (32) 

Here, an important remark is in order. The energy- 
domain formulation of Kirchhoff diffraction, Eq. (|3T))) . 
only assumes the validity of the Hclmholtz equation for 
a field in question, and is therefore applicable to a wide 
range of wave phenomena encountered, for instance, in 
acoustics, optics, and non-relativistic quantum mechan- 
ics. Equation (|32l) however relies on the particular re- 
lation, given by Eq. (0, between the energy-dependent 
Green function and time-dependent propagator, and is 
restricted to non-relativistic quantum mechanics only. 

Also, we note that the physical picture implied by 
Eq. ([32"T) is that of a particle traveling freely from q' 
to a point q! in the aperture, and then from q! to q. 
ft is therefore implicitly assumed that the aperture and 
points q and q' are chosen in such a way that every path 
q' — » q! > q' has no intersection with the screen other 
than at qi. In other words, Eq. (|32|) is only applicable 
to configurations in which apertures are not "shadowed" 
by S and are directly "visible" from the points q and q'. 



D. Diffraction in space and time 

We now note a striking similarity between (i) the 
Huygens-Fresnel expansion of the propagator for the 
problem of diffraction in time, Eq. (1181) , and (ii) the time- 
dependent formulation of Kirchhoff theory, Eq. (|32|). 
This leads us to a conjecture that both expressions are 
particular cases of a more general propagator expan- 
sion, describing quantum diffraction on apertures which 
themselves may vary in the course of time. Such time- 
dependent apertures are represented by a characteristic 
function x that depends on both the position qi along the 
dividing surface S and the instant t\ of the time interval 
(0,t). As before, x = x(qi!*i) i s allowed to take values 
between (perfect absorption) and 1 (perfect transmis- 
sion). In this case, we conjecture that the propagator, 
connecting two points q and q' on the opposite sides of 



the screen (such that s(q') < and s(q) > 0) in time t, 
is given by 

KM;t) = l I dh f dAfa *(«(<!!)) 
z Jo Jut 

( q - qi , qi - q' \ _ , s 

x tfo(q-qi;i-ti)x(qi;*i)^o(qi-q';*i)- (33) 

Equation ()33|) . which from now on we will refer as to 
Huygens-Fresnel-Kirchhoff (HFK) construction, consti- 
tutes the main result of the present paper. 

Our arguments supporting the conjecture, given by 
Eq. (|3"3")1 , are as follows. First, the propagator expansion 
is in accord with the Huygens-Fresnel principle. Second, 
the HFK construction correctly captures the physics of 
diffraction in time: In the special case that x is indepen- 
dent of qi and that the surface S is given by an (/ — 1)- 
dimensional plane, s(qi) = n-(qi — qo) with a unit vector 
n and some fixed vector q , Eq. (|3"3")) becomes equivalent 
to Eq. (fT8|) with x, x', and x\ replaced, respectively, by 
n • q, n • q', and n • qi. Third, in the case that x is inde- 
pendent of t\, Eq. p3p is nothing but a time-dependent 
formulation of Kirchhoff diffraction. 



III. DIFFRACTION OF SPATIALLY 
LOCALIZED WAVE PACKETS 

In this section, we apply the HFK construction, 
Eq. (|33[) . to wave functions that are initially localized in 
position space. We direct out attention to some example 
systems in one and two spatial dimensions. 

As before, we consider an absorbing screen that is de- 
fined by a surface S, Eq. (P2§1) . and a characteristic func- 
tion x- The latter, in general, is a function both of the 
position on S and of time. We further consider a quan- 
tum particle described at time t = by a wave function 
^(q'; 0), which is localized in a small neighborhood of a 
linear size a around a spatial point Q, i.e., ^(q'jO) ~ 
for |q' — Q| » ct. Here, we consider the case of a being 
small compared to the distance L between the point Q 
and the surface S. For concreteness, we take s(Q) < 0. 
Then, in accordance with Eqs. (QJ and (|33|) . the wave 
function $ at a point q, such that s(q) > 0, and time 
t > can be approximately written as 

*(q;t)~ I [ dh f d'qi^cu)) 

1 Jo JvJ 

/q-qi , qi-Q\ _ , , 

x U^ +- ^-J' Vs(qi) 

x iTo(q-qi;t-ti)x(qi;ii)*o(qi;*l), (34) 

where 

*o(qi;*i)=/ dV#o(qi-q';ti)*(q';0) (35) 
Jut 
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is the wave function of the corresponding free particle. 
Equation (|34p holds to the leading order in the small 
parameter a/L. 

A. One dimension 

As our first example we consider the Moshinksy prob- 
lem in one dimension (/ = 1) for a quantum particles 
initially described by the Gaussian wave packet 

(36) 

Here, Q and P represent, respectively, the average posi- 
tion and momentum of the particle, and a characterizes 
the wave packet dispersion in the position space. In our 
set-up, an absorbing shutter is placed at a point d, such 
that Q < d and a ^ L = d— Q. We are interested in the 
wave function ^>(x; t) at x > d and t > 0. 

A substitution of Eq. (|36[) . along with s(xi) = x\ — d 
and x = x(^i): mt o Eqs. (j3~4")l and ([33]) yields 

. . 1 /"* fx- d d-Q\ 

x K (x-d;t-tx) X {ti)%o{d;ti) (37) 

with 

Mx;t) = ( -^r^) 
\^Tt a ) 

and 

p ht 

Qt = Q + —t and lt = l + i ? . (39) 

m ma 

Evaluating the integral in the right-hand side of Eq. ([57)1 
numerically on obtains the wave function ^(x, t) at x > d 
and t > 0. 

Figure Q] shows the diffraction patterns calculated in 
accordance with Eq. (|37|) for different shutter protocols, 
x(ti). The initial wave packet, Eq. (|3"r!)) . is centered 
around Q — and has the average momentum P = 200 
and the position dispersion a = 0.1. Hereinafter, we use 
the atomic units, m = h = 1, in all numerical examples. 
The absorbing shutter is positioned at d = 8. Figure Q] 
shows the probability density ^(a;,^! 2 as a function of 
x for a fixed time, t = 0.05. Note that, for our choice 
of parameters, the position at the time t of the corre- 
sponding free classical particle is Qt = Q + Pt/m = 10. 
The figure shows four probability density distributions 
for (i) x(ii) = 1, representing the free-particle case, (ii) 
x(ii) = Q(ti — t ), corresponding to the perfectly ab- 
sorbing shutter first being in place and then suddenly 
removed at time to, (iii) x(^i) — ©(^o — t±), correspond- 
ing to the shutter closed instantaneously at to, and (iv) 
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FIG. 1: (Color online) The probability density |*(a;,t)| 2 , cal- 
culated in accordance with Eq. (|37p for various characteristic 
functions x(^i)i as a function of the position x. The initial 
wave packet is characterized by its position Q — 0, average 
momentum P — 200, and dispersion a = 0.1. The position 
of the shutter is given by d — 8, and the propagation time is 
t — 0.05. The shutter switching time to = 0.04 and the time 
window e = 5x 10~ 4 . All quantities are given in the atomic 
units, m = H = 1. 

x{ti) — 6(1— j*i — *q | /e) , representing a scenario in which 
the shutter is open only during a time interval of half- 
width e centered around to- Here we take to = 0.04 and 
e = 5 x 10~ 4 . Note that the position at time to of the 
corresponding free classical particle is Q + Pt /m = 8 
and coincides with the position of the shutter. 

Oscillations in the spatial dependence of the probabil- 
ity density, seen distinctly in Fig. [TJ are introduced by an 
instantaneous process of switching the shutter, and are, 
in fact, a manifestation of the diffraction-in-time phe- 
nomenon. It is well known that these oscillations become 
less pronounced and eventually disappear as one switches 
the shutter "continuously" over longer and longer time 
intervals [l(J [13 ■ Here we note that the HFK construc- 
tion provides a convenient framework for analyzing the 
disappearance (also known as apodization) of the diffrac- 
tion pattern for initial states that are localized in the 
position space. The phenomenon of the apodization of 
atomic beams, which correspond to initial states localized 
in the momentum space, has been previously addressed 
by different methods p7l[l8j. 

B. Two dimensions 

We now address diffraction of quantum wave packets 
in two dimensions on a (generally curved) shutter whose 
absorption properties may both depend on position and 
vary in time. We consider a quantum particle with the 
initial state given by 

*{x,y;0) =0i(x,O)0 2 (y,O), (40) 



where 



exp 



(C-Qj 

2oJ 



(41) 

for j = 1 and 2. Here, Q = (Qi,Q2) is the average po- 
sition and P = (P\,P2) average momentum of the par- 
ticle, and o i and a% characterize respectively the x- and 
y-component of the wave packet dispersion in the posi- 
tion space. As in the one-dimensional case, the initial 
distance between the particle and the shutter is assumed 
to be large compared to the spatial extent of the wave 
packet. 

Analogously to Eqs. (|38|) and (l39|) , the time evolution 
of a free particle is described by the wave function 



where 



*o(a,y;*) = <j>i(x,t) <fa{y,t) , 



(42) 



y t 



Q 



a < 





a > 



(a) 




CIO 



-8 



-7 



-6 



-1 



2 2 



exp 



-—t 

h 2m 



1 2 7j,t&j 



and 



p 

m 



and 



7i,t 



= l + i- 



(43) 



(44) 



for j = 1 and 2. Then, using this expression for \l/o in 
Eq. (|34[) and evaluating numerically the two-dimensional 
integral (over time and along a one-dimensional curve 
accommodating the shutter), we obtain the wave function 
^{x, y; t) at t > and at a point (x, y) on the side of the 
shutter opposite to Q. 

Figure [5] illustrates the effect of shutter curvature on 
the diffraction pattern. The set-up is schematically pre- 
sented in Fig.|2Ji. The initial state of the particle is given 
by a Gaussian wave packet, Eqs. (|4U|) and (|41[) . with the 
average position Q = (0,0), average momentum P = 
(200, 0), and dispersion at = a% = 0.1. The wave packet 
is incident upon a shutter, whose spatial geometry is de- 
fined by the curve s(xi, y±) — x\—d—ay\ = with d = 8. 
Depending on the sign of the parameter a, the shutter is 
concave (a < 0), flat (a — 0), or convex (a > 0). The 
time dependence of the shutter is given by the character- 
istic function x(Qii *i) = ©(1 — 1*1 — *o|/ e ) with t = 0.04 
and e — 5 x 10~ 4 . This corresponds to the case of a per- 
fectly absorbing shutter that is open only during a time 
interval of half-width e centered around to (cf. solid red 
curve in Fig.rjJ). The propagation time is set to t = 0.05, 
and the wave function is investigated inside a spatial re- 
gion defined by 8.9 < x < 11.1 and —1.6 < y < 1.6, 
shown schematically as a (light blue) rectangle in Fig. [5^. 
Note that the position of the corresponding classical par- 
ticle at time t equals Qt = Q + Pt/m = (10,0) and 



FIG. 2: (Color online) Diffraction in time on a curved shutter. 
Part (a) of the figure illustrates the set-up. The initial wave 
packet is characterized by its average position Q = (0,0), 
average momentum P = (200,0), and dispersion <j\ = 02 = 
0.1. The shutter is positioned along the curve s(xi, yi) — xi~ 
d — ayl — with d = 8, and its time dependence is specified 
by x(ti) = ©(I - 1*1 - to\/e) with to = 0.04 and e = 5 x 10" 4 . 
The propagation time is t = 0.05. The (light blue) rectangle, 
covering 8.9 < x < 11.1 and —1.6 < y < 1.6, outlines the 
spatial region, in which the probability density is computed. 
Parts (b), (c), and (d) of the figure show In \^(x, y; t)\ 2 for 
a = —0.3, 0, and 0.3 respectively. All quantities are given in 
the atomic units, m = ft = 1. 



coincides with the center of the rectangle. Then, the log- 
arithm of the probability density, In \^(x, y; i)| 2 , inside 
the rectangular region is shown for three different values 
of the parameter a: Fig. [2)3 corresponds to a = —0.3 
(concave shutter), Fig. [2t to a = (flat shutter), and 
Fig. [Hi to a = 0.3 (convex shutter). While the mean po- 
sition of the quantum particle coincides with that of the 
corresponding classical particle, the diffraction patterns 
are clearly different for the three choices of a. Compared 
to the diffraction pattern produced by the flat shutter, 
the concave shutter gives rise to a "divergent" pattern, 
while the convex shutter produces a "convergent" diffrac- 
tion pattern. 

Figure [3] addresses diffraction on a flat screen whose 
absorbing properties depend on time as well as vary in 
space. The set-up, schematically illustrated in Fig. 
As before, the initial wave packet is described by Q = 
(0,0), P = (200,0), a x = <J2 = 0.1, and the propaga- 
tion time is t = 0.05. Here we consider a flat screen, 
given by s{x\, y\) — x\ — d with d = 8, and study diffrac- 
tion patterns for three different characteristic functions 
X- As in Fig. [21 we present the logarithm of the prob- 
ability density, In \$?(x, y; t)\ 2 , inside a spatial region de- 
fined by 8.9 < x < 11.1 and —1.6 < y < 1.6, and illus- 
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FIG. 3: (Color online) Diffraction on a flat absorbing screen 
with a time-dependent aperture. Part (a) of the figure il- 
lustrates the set-up. Parameters of the initial wave packet 
are the same as in Fig. H Q = (0,0), P = (200,0), and 
<7i = (72 = 0.1. The propagation time is t = 0.05. The plane 
of the absorbing screen is defined by s(xi, yi) = x\ — d with 
d — 8. The (light blue) rectangle, covering 8.9 < x < 11.1 
and —1.6 < y < 1.6, outlines the spatial region, in which 
the probability density is computed. Parts (b), (c), and (d) 
of the figure show ln\^>(x,y;t)\ 2 for x(Qi>*i) equal, respec- 
tively, to 9(1 - \yi\/S), 6(1 - \ yi \/S)Q{l - \h - t \/e), and 
(1 - yl/5 2 - (ti - t f/e 2 ) with 5 = 0.05, to = 0.04, and 
e = 5 x 10~ 4 . All quantities are given in the atomic units, 
m = h — 1. 



trated by a (light blue) rectangle. Figure [3b shows the 
diffraction pattern produced by a slit in the absorbing 
screen and corresponds to x(qi;ii) = 0(1 — \Vi\/^) with 
S = 0.05. Figure [3J: represents a scenario, in which the 
slit of spatial half-width 5 is open only during a time in- 
terval of half-width e = 5 x 10~ 4 centered around time 
to = 0.04. The corresponding characteristic function is 
X(qi;*i) = ©(! - M/Wl - 1*1 - t Q \/e). Figure^ 
shows the diffraction pattern obtained for the character- 
istic function x(qi;ii) = 6(l - yl/5 2 - (t x - t ) 2 /e 2 ). 
In this case, the width of the slit continuously increases 
from to 25 during the time interval to — e < t\ < to, 
and then shrinks back to during to < t% < to + e. 

It is interesting to observe some similarity of, on 
one hand, the diffraction patterns of Figs. [3J: and [5J1, 
and, on the other hand, patterns of Fraunhofer diffrac- 
tion of light on, respectively, rectangular and elliptic 
apertures in opaque two-dimensional screens (e.g., see 
Figs. 8.10 and 8.12 in Ref. @). The similarity becomes 
even stronger when one takes into account the space- 
time, (2/i,ti) representation of the apertures leading to 
the diffraction patterns shown in Figs. [3J: and [3Ji. In- 
deed, the aperture defined by the characteristic function 



X(qi;*i) = ©(I - M/<*)6(1 " 1*1 " to\/e) can be rep- 
resented as a rectangular opening in the (yi,t\) coordi- 
nates, while the aperture corresponding to x(qi;ti) = 
0(l — yl/5 2 — (ii — to) 2 /e 2 ) has the shape of an ellipse in 
the same coordinates. This observation suggests that /- 
dimensional quantum dynamics in the presence of time- 
dependent diffracting obstacles may be used for modeling 
diffraction in (/ + l)-dimensional systems with station- 
ary, time- independent obstacles. 



IV. CONCLUSIONS 

In this paper we have revisited the phenomenon of 
diffraction in non-relativistic quantum mechanics. More 
specifically, we have considered the problem of a quantum 
particle passing through an opening (or a set of openings) 
in a perfectly absorbing screen. Having assumed the 
validity of the Huygens-Fresnel principle, we have con- 
structed an expansion of the quantum propagator that 
connects two spatial points lying on the opposite sides of 
the screen. The expansion represents the full propagator 
as a sum of products of two free particle propagators, one 
connecting the initial point and the other connecting the 
final point to a point in the opening. The construction 
holds in the cases of curved (convex or concave) screens 
and for openings whose shape changes in time. Conse- 
quently, our approach allows one to analyze the quantum 
phenomena of diffraction in space and diffraction in time, 
as well as the interplay between the two. 

In order to illustrate the method, we have used our 
propagator expansion to calculate diffraction patterns for 
an initially localized wave packet passing through vari- 
ous spatio-temporal diffraction screens in one and two di- 
mensions. Thus, in two dimensions, we have investigated 
the effect of screen curvature on the diffraction pattern 
by analyzing diffraction in time produced by absorbing 
parabolic shutters and demonstrating how a convex (con- 
cave) shutter gives rise to effective focusing (defocusing) 
of the wave function. We have also studied diffraction of 
spatially localized wave packets passing through holes of 
time-dependent size. We have shown that the shape of 
a diffracted wave function is largely determined by the 
space-time geometry of the hole, suggesting the use of 
time-dependent obstacles in /-dimensional systems for 
modeling quantum diffraction on stationary obstacles in 
(/ + l)-dimensional systems. 

The approach, developed in this paper, is applicable 
to quantum diffraction on perfectly absorbing screens, 
and, therefore, complements the method of Brukner and 
Zollinger [l5[ valid for perfectly reflecting screens. In 
laboratory experiments, however, diffraction screens are 
typically neither perfectly absorbing nor perfectly reflect- 
ing, and the two methods have to be used in combination. 
It would interesting to develop such a unified approach 
and to test its predictions against results of experimental 
measurements. 
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Appendix A: Derivation of Eqs. (|10|) and 



where W stands for the Whittaker function. The latter 
can be written as (see page 341 in Ref. [25[, or section 
13.18(h) of Ref. 



W-i.(( 2 ) = exp(C 2 /2)erfc(C) 



(A6) 



with "erfc" denoting the complementary error function. 
Substituting Eqs. (|X5|) and ([M)) into ([M)) we obtain 



In one dimension the free-particle Green function, 
Go(/; E) with I > 0, is given by 

/>oo 

G {l;E)=L[K ](l;s) = / dte- st K (l;t) 
Jo 



in 



dt 



2mh J y/t 



exp 



ml 
2AM, 



st 



to / /2ms 



(Al) 



Here, s = E/(ih) is implied, cf. Eq. ((5|), and the reader 
is referred to the formula 3.471.15 of Ref. [24[ for the last 
integral. 

We now consider the integral 

X= [ dtiKo(x-xi;t-ti)Ko(xi-x';tx), (A2) 
Jo 

where x, x' , and x\ satisfy Eq. ([7]). The Laplace trans- 
form of X is given by 

C[X] = C[Kq](x — x\; s) L[Kq](xi — x';s) 
in 



2ihs 



exp \ ~(x — x') 



2ms 
ih 



rn 
2ihs 



£[K }(x-x'; S ). 



(A3) 



Then, using the fact that 1/y/s — £[l/\/irt] and perform- 
ing the inverse Laplace transform, we obtain 



1 = 



dt 1 



2-1X1% J y/t — t\ 

* dti 



K (x - x';tx) 

m(x — x') 2 



exp 



2ihh 



(A4) 



The integral in the last line of Eq. (|A4[) is given by the 
formula 3.471.2 of Ref. [11, 



dti 



exp , 

y/hit-h) V *1 



2t i'i V t 



(A5) 



(A7) 



1- = tt^t erfc ( \ / — — (x — x') 
2ih \\ 2ikV ' 



A comparison of the right-hand sides of Eqs. (|A2[) and 
(IA7|) completes the proof of the identity (|9|) with the 
function u chosen in accordance with Eq. (jlip . 



In order to prove that Eq. (1101) offers an alternative 
choice for the function u we consider 



(A8) 



where X is given by Eq. (|A2[) . and r\\ and 772 are two 
arbitrary, generally complex numbers. On one hand, X' 
can be evaluated by a direct substitution of Eq. (|A2[) into 
Eq. ([Ml) which leads to 



. 7777 f ( X — X\ X\ — X 1 



x K Q (x - xi;t - ti)K (xi - x';ti) . (A9) 



On the other hand, using Eqs. (|M|) . (TM)) . and (fATj) . we 

have 



d d 
^xl + ^jrj ) C[K ]{x-x';s) 



2ihs V dx dx' 



'-(m - m)£[ K o](x - x';s) , 



(A10) 



and therefore 



X' = —( m - V2 )K (x-x';t). (All) 

h 



Finally, comparing Eqs. (|A9I) and (IA11[) . and introducing 
a complex number r\ = 7/1/(7/1 — 7/2)1 we arrive at the 
identity ^ with the function u given by Eq. (ITU1) . 
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